Fermionic Superfluid from a Bilayer Band Insulator in an Optical Lattice 
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We propose a model to realize a fermionic superfluid state in an optical lattice circumventing the 
cooling problem. Our proposal exploits the idea of tuning the interaction in a characteristically low 
entropy state, a band-insulator in an optical bilayer system, to obtain a superfluid. By performing 
a detailed analysis of the model including fluctuations and augmented by a variational quantum 
Monte Carlo calculations of the ground state, we show that the superfluid state obtained has high 
transition temperature of the order of the hopping energy. Our system is designed to suppress other 
competing orders such as a charge density wave. We suggest a laboratory realization of this model 
via an orthogonally shaken optical lattice bilayer. 
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Quantum emulation of interesting condensed matter 
Hamiltonians using ultra cold atom systems holds much 
promise. [1-5] The pace of experimental progress has 
been impeded by key problems which include simula- 
tion of electromagnetic (gauge) fields, removal of entropy, 
etc. While the former has seen a spectacular recent 
progress [6-10], the long standing "cooling problem" of 
trapped lattice fermions has been more difficult. [11] 

The cooling problem has been addressed in various 
ways. One approach has been to find schemes to 
"squeeze out" the entropy. [12, 13] Others include exploit- 
ing metastable states[14], using properties of the states 
(such as the Neel state) to develop cooling protocols[15] 
(see ref. [11] for a review). A recent notable proposal is 
to use an additional beam that helps to enlarge the re- 
gion where a desired state is stabilized. [16] Despite this, 
to the best of our knowledge, an interesting many body 
state such as an anti-ferromagnet is yet to be realized 
in an optical lattice, while some signatures of fermionic 
superfluidity have been reported. [17] 

A different strategy would be to create a characteristi- 
cally low entropy state in a large region of the trap, and 
to tune a parameter (such as the interaction) that drives 
this low entropy region into an interesting many-body 
state. The desiderata of such a many-body state include 
(i) high characteristic temperature scale (ii) stability over 
other "uninteresting" competing states. Here we suggest 
the use of a band insulator, a characteristically low en- 
tropy state in which we tune an attractive interaction 
to produce a fermionic superfluid. Band-insulator super- 
fluid transitions have been investigated earlier in other 
contexts. [18, 19] In fermionic cold atom systems, moti- 
vated by experimental work cited above[17], superfluid- 
band insulator transition engendered by increasing the 
lattice depth with concomitant multiband effects have 
been discussed. [20-24] In contrast to these works, our 
proposal aims to obtain a superfluid in a deep lattice. 

We propose and study a bilayer band insulator that 
undergoes a transition to a superfluid upon tuning an at- 
tractive interaction, attaining the above desiderata. The 




Figure 1. Dispersion of the bilayer band insulator - 

Bilayer dispersion before (Left) and after (Right) interlayer 
hybridization. Inset show schematic density of states. 



model is designed so that competing phases, such as 
charge density wave (CDW), are avoided. This is demon- 
strated by a detailed analysis including Gaussian fluc- 
tuations, and variational Monte Carlo simulations. We 
show that a "high-temperature" superconducting phase 
is possible in this system by estimating the Berezinski- 
Kosterlitz-Thouless transition temperature Tbkt- In 
a regime of parameters, the system shows interesting 
physics such as "pseudogap phenomenon", even at high 
temperatures. We suggest a possible route to realize this 
in an optical lattice. 

Our proposal for the realization of a two dimensional 
spin- 1 fermionic superfluid state hinges on a bilayer band 
insulator. The configuration consists of two layers A and 
B, both of which have the same lattice structure (such as 
a square or triangular lattice) and a 2D Brilliouin zone 
(see Fig. 1). The crucial ingredient is that the in-plane 
energy dispersion in the two layers are of opposite sign, 
i. e., SA(k) = — Eg(k) = e(k) for all k in the Brilliouin 
zone. Inter-layer hopping is described by a hybridization 
function h(k) which is such that h(k) and e(k) never van- 
ish simultaneously. The kinetic energy of the system is 

thus described by U K = J2 k a £ ( k ) ( a L a *^ - &L&fc<r) + 

(^(tya^bka + h. c.V where a-s and b-s are spin-| 
fermion operators corresponding to A and B layers re- 
spectively (cr =t,l is the spin). This leads to two (spin 
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degenerate) bands 

n K = ^2e(k)(c{ r7 c ktT -d{ (T d k<T ), (1) 

ka 

where c and d are respectively "conduction" and 
"valance" band fermion operators, with e(fe) = 
\J 'e(fe) 2 + \h(k)\ 2 . With a fermion density of one par- 
ticle per site on both A and B layers, the ground state is 
the filled valance band, i. e., a band insulator. We now 
introduce a local attractive interaction with strength U 
(which may be tuned by a Feshbach resonance [2 5]) as 

n U = ~U ( a It a 4 a 4 a *T + &l t &ljA4&;t) • (2) 
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We show that a superfluid can be generated by start- 
ing from the band insulator (U ~ 0) and adiabatically 
increasing the magnitude of the attractive interaction. 
To be specific, we choose a particular model to illus- 
trate the idea (see below for a possible laboratory re- 
alization of this model). The bilayer system has both 
A and B layers which are square lattices (of unit lattice 
spacing) with nearest (t) and next nearest (t') hopping 
(i' = t/10 throughout) such that SA(k) = — £s(fe) = 
e(k) = — 2i(cos k x + cos k y ) — At' cos k x cos k y . The hy- 
bridization function h(k) = —th captures the hopping 
from adjacent A and B sites. The resulting band struc- 
ture is of the form in eqn. (1) and has an energy gap 
e g = 2th- The band insulator obtained with one par- 
ticle per site per layer has a Cooper instability [26] at a 
non-zero critical value U c of U unlike in a metal where 
U c is zero. We find that i « ^E ke Jk) ~ I ln ^ 
(see Fig. 2), where N is the number of sites per layer. 
The instability owes to the fact that at sufficiently large 
U, it becomes feasible for a pair of fermions of opposite 
spin to be promoted to conduction band where they can 
"sample" the attractive interaction (see Fig. 2(a)), even- 
tually forming a bound state. We find that for U > U c , 
the binding energy of the pair goes as (U — U c ), which 
may be contrasted with exponentially small value usu- 
ally found[26] in a system with a Fermi surface. The 
physics of such strong binding owes to the modification 
of the density of states at the band edges engendered by 
the hybridization. The resulting joint density of states of 
particle-hole excitations is strongly enhanced (see Fig. 1), 
g(e) ~ — e g , and it is this large enhancement that 

provides for the strong binding as in other contexts. [27] 
Consequently we expect the system to also possess high 
transition temperatures making it attractive for experi- 
mental realization of an optical lattice superfluid. 

We now study the properties of the lattice superfluid 
state using functional integral techniques. [28, 29], by in- 
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Figure 2. (a) Schematic of Cooper instability, (b) The critical 
value U c that induces Cooper instability. 

troducing the action 

k,oc 

k q,a 

(3) 

where — Gq l {k) = (—ik n + ae(k) — fi), fi is the chem- 
ical potential, k — (ik n ,k), (iqt,q) with ik n (iqi) — 
(2n + 1 )vr/(3 (2£w/fi) being Fermi(Bose) Matsubara fre- 
quencies, and |3 = 1/T is the inverse temperature. 
We have introduced Grassmann numbers "0 QCT where the 
flavour label a = ±1 stands, respectively, for A and B 
layers, and P*(q) = £ fe $^(g + k)ip*^(-k). 

Possibility of a superfluid state is investigated by in- 
troducing Hubbard-Stratanovich pair fields A a (g) to de- 
couple the interaction term[30] in eqn. (3). The fermions 
are then integrated out to obtain an action 5 [A] solely 
for A a (q). The uniform saddle point A„(g) = A6 q fi, 
where A is the superfluid order parameter, gives the gap 
equation 

1 1 tanhM 

U N ^ 2E(k) ' 1 ' 

where E{k) = J e(k) 2 + A 2 . The model we consider 
here has particle hole symmetry which forces fj, = when 
the occupancy is one fermion per site; thus a separate 
number equation to determine /i is obviated. Fig. 3(a) 
shows the evolution of zero temperature ground state 
with increasing U. At U c , a quantum phase transition 
occurs ushering in a superfluid state where A behaves 
as y/U — U c and monotonically increases with increasing 
U. The superfluidity arises from promotion of fermions 
to the conduction band; indeed, n = -A* J2k( c kcr c ko-) ( see 
eqn. (1)), the number of fermions promoted to the con- 
duction band increases from zero at U c as n ~ (U — U c ). 

Having established the superfluid ground state and its 
physical underpinnings, the natural question is regarding 
the magnitude of the transition temperature of the su- 
perfluid obtained. Two effects that destroy superfluidity 
are pair breaking and phase fluctuations. Indeed, long 
wavelength phase fluctuations render our 2D superfluid 
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Figure 3. Dependence of zero temperature (a) superfluid or- 
der parameter A and (b) "carrier density" n on U . 




Figure 4. Evolution of (a) superfluid density p s and (b) 
Leggett mode gap ljl as a function of U at zero tempera- 
ture. 



lacking in true long range order at finite temperatures, 
i. e., in the Kosterliz-Thouless phase. The temperature 
scale of pair-breaking T& is set by the lowest temperature 
at which the saddle point value of A vanishes, and can 
be obtained by solving for the temperature in the gap 
equation (eqn. (4)) with A — 0. To investigate the role 
of phase fluctuations below T& and estimate the transi- 
tion temperature, we study the fluctuations at the Gaus- 
sian level [30] by expanding the action S [A] about the 
saddle point with A a (q) = A%(q) + A(( a (q) + i6 a (q)) 
where C, a {q) and a {q) are real fields that represent, re- 
spectively, the amplitude and phase fluctuations in each 
layer. The fluctuations in each layer are coupled; a more 
natural "normal mode" description is in terms of sym- 
metric and anti-symmetric linear combinations of these 
modes. For example, there are two phase modes - the 
symmetric mode 9 s {q) ~ (#+(<?) + 9-(q)), and the anti- 
symmetric mode a (q) ~ (#+(?) — 0_ (<?)), and there are 
two amplitude modes with similar definition. 

We find that both the amplitude modes are gapped, 
while the symmetric phase mode is gapless and the anti- 
symmetric phase mode is gapped. Interestingly, the 
gapped anti-symmetric phase mode is analogous to the 
Leggett mode in multi-band superconductors. [31, 32] We 
obtain the following effective action for the phase modes 
by integrating out the amplitude modes S[9 S 
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u)i,0£j , where r is the imaginary time, r 

is the position on the 2D plane, k-s, p-s and u>l are de- 
termined by the saddle point solution. The most impor- 
tant parameters [30] in this action are the phase stiffness 
of the symmetric phase mode p s , which is the superfluid 




Figure 5. (a) Determination of Tbkt from the temperature 
dependence of p a . (b) Dependence of Ta and Tbkt on U. 
Note the high transition temperature, and the large pseudo- 
gap regime. Here tf, = 0.5t. 



density, and ojl, the gap (or mass) associated with the 
anti-symmetric Leggett mode. Fig. 4(a) shows the depen- 
dence of the zero-temperature superfluid density p s on U. 
For U > U c , p s ~ U —U c , and has the same behaviour as 
the number of fermions excited to the conduction band. 
With increase of U, p s attains a maximum, and suffers a 

.2 

fall at larger values of U. For U » i, we find that p s ~ jj. 
This owes to the fact that the system undergoes a "BCS- 
BEC" crossover with increasing U, and p s is determined 
by the hopping amplitude of the bosonic fermion pair at 
large U which is ~ t 2 /U. The variation of with U 
is shown in Fig. 4(b), the key point to be noted is that 
in the regime where p s is largest, the Leggett mode has 
a large gap, and does not participate in the low energy 
physics. 

The discussion above allows the estimation of the 
Kosterlitz-Thouless transition temperature Tbkt- We 
obtain p s as a function of T via our functional formu- 
lation using the saddle point value of A. Using the 
relationship [33] that p s {Tbkt) — 2Tbkt , we arrive at 



the transition temperature (as shown in Fig. 5(a)) plot- 
ted in Fig. 5(b), which also shows the temperature Ta 
associated with pair breaking obtained from eqn. (4). We 
see that the maximum value of Tbkt is of the order of 
the lattice hopping amplitude t, and in this sense we ob- 
tain high temperature superfluidity in the regime where 
the crossover to the BEC state takes place. The BCS 
side (U > U c ) is also a robust superfluid owing to the en- 
hancement obtained by the divergent density of states. 
The transition temperature of the system can also be af- 
fected by the vortex core energies [34] , but these effects 
will be unimportant in this system. Another attractive 
aspect of this system is that one expects to see large pseu- 
dogap features even at high temperatures (see Fig. 5(b)), 
and thus interesting physics can be investigated in opti- 
cal lattices even if the average entropy of the system is 
not small. 

The effect of quantum fluctuations are likely impor- 
tant due to the reduced dimensionality. [35] To ensure 
that quantum fluctuations only have a quantitative role, 
and to ensure that there are no competing orders such 
as a CDW intervening, we conducted a detailed varia- 



4 




02468 02468 



U/t U/t 

Figure 6. Results of variational Monte Carlo calculations. 
Dependence of (a) variational parameters g and As and (b) 
superfluid order parameter $, on U . th = 0.5t. 




Figure 7. Orthogonally shaken bilayer optical lattice. The top 
layer A is shaken in the x direction, while the bottom layer B 
is shaken in the y direction. By an appropriate choice of laser 
intensities, amplitude and frequency of the shakes, the band 
dispersion in the layers can be made to have opposite signs to 
each other. The layer hybridization can be controlled by the 
distance between the layers. 

tional Monte-Carlo calculation of the ground state. [30] 
Our variational ground state |^) = g D \As, Acdw) bcs 
is constructed by introducing both the superfluid pair or- 
der As, and a commensurate (w, it) charge density wave 
order parameter Acdw, and obtaining the BCS state 
\Ag, Acdw) bcs- The Gutzwiller parameter g (> 1) 
that promotes double occupancy (D is the operator that 
counts the number of doubly occupied sites) introduces 
quantum fluctuations of the local phase. Two key results 
of our detailed study are: (i) For all values of U within 
the range considered here, the optimal value of Acdw is 
zero, i. e., there is no competing order that intervenes and 
hence the superfluid state is stable, (ii) quantum fluctua- 
tions do not change the qualitative aspects of the results. 
Indeed, for the parameter values shown in Fig. 6(a), we 
find U c ~ 3.2 is expectedly larger than the value of 2.3 
from the saddle point analysis. The variational param- 
eter As (Fig. 6(a)) and the superfluid order parameter 
$ (Fig. 6(b)), which measures the amplitude of injecting 
a pair at a large distance away from the point of its re- 
moval, has precisely the behaviour as expected from the 
saddle point analysis. 

Our proposed scheme can be realized by an "orthogo- 
nally shaken bilayer" depicted in Fig. 7. It was argued 
in [36] that on introducing a shake of the optical lat- 
tice, the amplitude and sign of the hopping can be con- 
trolled. It was shown that if K (an energy scale) and 
v are respectively the amplitude and frequency of the 
shake, the effective hopping amplitude t e ff — tJo(v) 



where Jo is the Bessel function and t is the hopping in 
the absence of the shake. This phenomenon has not only 
been observed experimentally [37], but has been recently 
used to study many interesting quantum phases [38, 39] 
with further proposals for the generation of topological 
insulators. [40, 41] 

Our proposed experimental set up consists of two adja- 
cent optical square lattices. The top layer A is obtained 
by interfering two sets of counter propagating laser beams 
in the x and y directions; the x and y beams are non- 
interfering. The relative phase of the two x laser beams 
is modulated so as to obtain a shake, and intensity of the 
x laser beams and amplitude of the modulation can be 
chosen such that — 1£ — ty — t, i. e., the hopping along 
the x direction has an opposite sign to that in the y direc- 
tion. In the layer B, the beams along the y-direction are 
shaken so that t£ = — ty = t. This provides a realization 
of a system with e^{k) = — £b(&)- The hybridization 
of the two layers can be controlled by the distance be- 
tween the two layers. This can be achieved by using 
vertically confining beams as in ref. [42], and creating 
the two layers by "optical copying" . Optical copying will 
entail splitting the x and y laser beams of the A layer 
and focusing the split beams just below the A layer to 
produce the B layer such that the x beam of the A layer 
plays the role of the y beam of the B layer, and so on. 
This laboratory realization of our proposal may require 
optics techniques that have been used in the making of 
quantum gas microscope. [43] The confining trap poten- 
tial is to be designed such that a large region near the 
trap center will be in a band insulating state, with the ex- 
cess entropy trapped in regions at the periphery. Tuning 
of the attractive interactions should now drive the cen- 
tral band insulating region to the superfluid state. We 
hope that this work stimulates experimental research on 
realizing such a bilayer band insulator system, even by 
routes other than our proposal. 
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